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Abstract 

Determining the existence of interconnected responses among life-history traits 
and identifying underlying environmental drivers are recognized as key goals 
for understanding the basis of phenotypic variability. We studied potentially 
interconnected responses among senescence, fecundity, embryos size, weight of 
brooding females, size at maturity and sex ratio in a semiterrestrial amphipod 
affected by macroscale gradients in beach morphodynamics and salinity. To this 
end, multiple modelling processes based on generalized additive mixed models 
were used to deal with the spatio-temporal structure of the data obtained at 10 
beaches during 22 months. Salinity was the only nexus among life-history traits, 
suggesting that this physiological stressor influences the energy balance of 
organisms. Different salinity scenarios determined shifts in the weight of brood- 
ing females and size at maturity, having consequences in the number and size 
of embryos which in turn affected sex determination and sex ratio at the popu- 
lation level. Our work highlights the importance of analysing field data to find 
the variables and potential mechanisms that define concerted responses among 
traits, therefore defining life-history strategies. 
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Introduction 

Environment as an ecological actor causes variations in 
phenotype (phenotypic variability) and therefore in the 
performance of individuals (Pigliucci 2001). The physio- 
logical interdependence among life-history traits con- 
strains and adjusts plastic responses, determining that the 
effects of a certain environmental pressure on one trait, 
has an effect on the others (Ricklefs and Wikelski 2002). 
Therefore, elucidating which variables are responsible for 
the observed phenotypic variability and the linkages 
among traits has been recognized as a key goal for under- 
standing the basis of phenotypic variability in natural 
populations (Ricklefs and Wikelski 2002; Isaksson et al. 
2011). 

Widely distributed species that are subject to a broad 
range of conditions represent useful frameworks for 



evaluating environmental effects on phenotypic variability 
(Joshi et al. 2001; Defeo and McLachlan 2005). Water-land 
transitions in the evolutionary history of the talitrid amphi- 
pods (Hurley 1968; Spicer et al. 1987), the variety of habi- 
tats in which they live and the large geographical 
distribution that they exhibit (Ramus and Forward 2012), 
made them key cases for this kind of study. Particularly, 
those that inhabit sandy beaches evolved to persist in this 
harsh environment, where life histories are considered to be 
mainly shaped by physical factors (Defeo and McLachlan 
2005). In this context, a morphodynamic spectrum from 
narrow and steep beaches (reflective beaches) to wide and 
flat ones (dissipative beaches) (Short 1999) affects sandy 
beach macrofauna, including amphipods (Defeo and 
Gomez 2005; McLachlan and Brown 2006). In addition, tal- 
itrid amphipods are affected by desiccation and osmotic 
constrains (Truchot 1990; Morritt and Spicer 1998), as they 
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are exposed to semiterrestrial conditions and salinity fluctu- 
ations (Gomez and Defeo 2012; Ramus and Forward 2012). 

Uruguayan sandy beaches cover a wide morphodynam- 
ic range along the Rio de la Plata estuary (the widest 
worldwide) and the Atlantic Ocean, providing an ideal 
scenario to study how the environment affects life-history 
traits. This is particularly true in the talitrid amphipod 
Altantorchestoidea brasiliensis (Fig. 1), which inhabits the 
whole morphodynamic spectrum and, at the same time, 
is affected by the estuarine gradient generated by the Rio 
de la Plata. We used field data of this talitrid amphipod 
and supervised machine learning methods with mixed 
models to elucidate the environmental drivers that affect 
the response of multiple life-history traits and to under- 
stand the involved mechanisms. 

Material and Methods 

Study area and field sampling 

The Uruguayan coastline stretches over 320 km along the 
north-eastern bank of the Rio de la Plata estuary and 
210 km of adjacent Atlantic coastline (Lercari and Defeo 
2006). Ten sandy beaches were sampled along this 
temperate microtidal coast, covering a wide range of 
morphodynamic states along the estuarine gradient. 
Environmental and biological samples were obtained 
every 2 months from July 1999 to April 2001. Three tran- 
sects were set up perpendicular to the shoreline and 
spaced 8 m apart. On each transect, sampling stations 
were carried out every 4 m from the base of the dunes to 
the seaward limit of the A. brasiliensis distribution, until 
at least two consecutive stations yielded no amphipods. 
At each station, a sheet metal cylinder (27 cm diameter) 




Figure 1. The talitrid amphipod Altantorchestoidea brasiliensis. 
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was used to remove the sediment up to a depth of 
40 cm. Samples were field-sieved through a 0.5 mm 
mesh; samples with coarse sand were carefully sorted in 
the laboratory to retain all organisms. At each sampling 
station, sediment samples were collected with a corer 
(6.5 cm diameter) for determination of mean grain size, 
sand moisture and organic matter content. Measurements 
of sediment temperature and sand compaction were 
made, in the latter using a piston pocket penetrometer 
(Herrick and Jones 2002). Salinity and water temperature 
were measured. Wave height and wave period were deter- 
mined visually using a stopwatch. Beach slope was esti- 
mated by Emery's profiling technique (Emery 1961). 
Beach width was measured as the distance between the 
base of the dunes and the lower swash level. 

In the laboratory, collected amphipods were measured 
from the tip of the cephalon to the end of the telson, 
using binocular microscopes equipped with micrometrical 
ocular lens, calibrated with objective micrometres. Indi- 
viduals were sexed based on the presence of copulatory 
appendages in males, and the presence of oostegites, with 
or without setae, in females (Marques et al. 2003). 
Females with setae in the oostegites were considered as 
mature, as they had, have or will have embryos. We con- 
sider as brooding females only those carrying embryos. 
Individuals with no secondary sexual dimorphic features 
were classified as juveniles. Embryos were separated from 
the brood pouch, counted and measured. Afterwards, all 
individuals were dried at 50°C for 24 h and weighted 
to the nearest 0.0001 g. We estimated abundance (indi- 
viduals per strip transect; ind-m~') (Brazeiro and Defeo 
1996) discriminated by size class for males, females and 
the whole population (adults and juveniles). Eleven 
size classes of 1.00 mm (from >2.50 to <13.50 mm) were 
considered. 

Modelling process 

We modelled abundance by size class, number of embryos 
per female, brooding females dry weight, embryos size, 
proportion of mature females and proportion of females, 
in relation to environmental variables, using generalized 
additive mixed models (GAMM) (see details in Table 1). 
To do so, the gamm4 R package was used (Wood 2011). 
A categorical site variable and the date of each sampling 
event were included as random intercepts in all models, 
to adjust the variability among sites and times without 
consuming a great amount of degrees of freedom. All 
potential predictors were included without theoretical 
constrains, being the parsimony of the models (evaluated 
using Akaike's Information Criterion), the only basis to 
retain them. Variables were considered as smooth terms 
using penalized regression splines with degrees of freedom 
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Table 1. Details of Generalized Additive Mixed Models used for mod- 
elling life-history traits of Altantorchestoidea brasiliensis. In each 
model, all potential predictors were included without theoretical con- 
strains; being the parsimony of the models, the only basis to retain 
them. 
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up to 3 (Hastie and Tibshirani 1990; Barboza et al. 2012; 
Gomez and Defeo 2012). Residuals were checked for all 
models. 

Results and Discussion 

Abundance models for males and females retained size 
class, and salinity as the best explanatory predictors of the 
observed variability (Fig. 2A-D). The partial response of 
log- abundance in relation to size class (considered as a 



proxy of age, see Gomez and Defeo 1999) showed that 
senescence (as a constant decay) in A. brasiliensis starts at 
ca. 9 mm both for males and females (Fig. 2A, C). In the 
particular case of females, this result coincides with the 
size at which 50% are in a brooding condition according 
to the observations of Gomez and Defeo 1999 and the 
results for female's maturity obtained herein (Fig. 4A). 
This suggests that A. brasiliensis females might be repro- 
ducing only once in this temperate region, increasing its 
mortality as a consequence of the large physiological 
energy invested in reproduction ("reproductive senes- 
cence", see Kirkwood and Austad 2000). Salinity was the 
single and main environmental predictor of adult abun- 
dance. Both sexes exhibited their highest abundance at 
salinities near 15, consistently with laboratory and field 
observations in other supralittoral amphipods (Steele and 
Steele 1991). The final model of total abundance (juve- 
niles and adults) retained beach slope in addition to size 
class and salinity (Fig. 2E-G) as explanatory variables. 
The retention of beach slope only in the total abundance 
model remarks the effects of beach morphodynamics on 
juveniles. The steep slopes in microtidal reflective beaches 
prevent unpredictable water intrusions in the supralittoral 
zone due to storms and winds, protecting juveniles that 
are incapable of escaping towards the dunes (Gomez and 
Defeo 1999; Defeo and Gomez 2005). 

The females fecundity model only retained individual 
weight as the main predictor (Fig. 3A), indicating that 
environmental conditions do not have a direct effect on 
the number of embryos produced. The increase in fecun- 
dity with individual weight reveals the importance of con- 
sidering body weight as proxy of the available energy to 
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Figure 2. Gaussian GAMM of abundance 
(log-transformed) by size class: females (A, B), 
males (C, D) and total population (adults and 
juveniles) (E, F, G). Partial effects of retained 
predictors (presented in order of importance) 
are shown. The numbers in they-axis labels 
indicate the estimated degrees of freedom of 
the smooth curve. Marks on the x-axis indicate 
the measured values of each variable. 
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Figure 3. Poisson GAMM for the number of embryos per female (A), 
and Gaussian GAMM for brooding females body dry weight (B, C) 
and for embryos size (D, E), showing partial effects of retained 
environmental predictors. The numbers in the y-axis labels indicate 
the estimated degrees of freedom of the smooth curve. Marks on the 
x-axis indicate the measured values of each variable. 



invest in reproduction. In addition, when the response 
variable was body weight of brooding females, individual 
size and salinity were the relevant variables (Fig. 3B, C), 
showing once again the importance of salinity in shaping 
the life history of A. brasiliensis. Species' preference for 
salinities close to 15-20 (Gomez and Defeo 2012) suggests 
that the decrease in dry weight towards fully freshwater 
conditions may be due to the energy expenditure required 
to allow an adequate osmoregulation of females. Metabolic 
rate may increase as a result of osmotic stress inflicted by 
the decline in salinity, decreasing the efficiency to absorb 
nutrients and generate reserves (Normant and Lamprecht 
2006). This becomes clear if we consider that ionic regula- 
tion is a steady state that requires large amounts of energy 
to be accomplished (Schmidt-Nielsen 1997). Furthermore, 
GAMM results revealed a shift in the size at which brood- 
ing females reach the same weight at low salinities, in 
comparison with those that inhabit salinities >15. Thus, 
for a given size, brooding females exhibit lower weights in 
low-salinity environments (Fig. 3C and Fig. 6). 

Embryos size decreased with salinity (Fig. 3D). This, in 
conjunction with the fact that females have a limited 
space in their pouches, allow us to consider the existence 
of a trade-off between embryos size and fecundity, taking 
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Figure 4. Binomial GAMM (logit scale) of maturity, showing partial 
effects of retained predictors on mature/nonmature females: (A) size 
and (B) salinity. The numbers in the y-axis labels indicate the 
estimated degrees of freedom of the smooth curve. Marks on the x- 
axis indicate the measured values of each variable. 



into account that heavier females have more embryos in 
environments with intermediate and high salinities. Thus, 
two reproductive strategies could exist: (1) low number 
of large embryos at low salinities; and (2) high number of 
small embryos at intermediate/high salinities (Fig. 6). One 
potential explanation is that females that occur in environ- 
ments where they have an efficient osmoregulation (salini- 
ties >15) exhibit higher weight and energy availability, 
being able to produce more embryos. By contrast, females 
in environments where they have a large expenditure in 
osmoregulation (salinities <15) invest their energy in size 
but not in number of embryos. The production of larger 
embryos in low-salinity systems allows future juveniles to 
invest more energy in osmotic control once they become 
independent from the parental care. An alternative expla- 
nation is the osmotic water uptake that embryos could 
present at low-salinity waters (Moran and McAlister 
2009). Moreover, embryos size decreased with temperature 
(Fig. 3E). Gammarus duebeni and G. insensibilis also pro- 
duce bigger embryos during cold months and smaller ones 
in the summer (Sheader 1983, 1996). The physiological 
bases for these observations are not understood, but tem- 
perature has a well-documented negative correlation with 
egg size in many marine taxa, both within and among 
species (Levitan 2000; Moran and McAlister 2009). 

Size at maturity in females (Fig. 4A) was affected only 
by salinity (Fig. 4B). This life-history trait exhibited a 
shift towards salinities <20. Thus, 50% of females reached 
maturity at larger sizes or, for a given size, the proportion 
of mature females decreased in low-salinity environments 
(Fig. 6). This coincides with the shift caused by salinity in 
the size at which females reach a certain weight, suggest- 
ing a metabolic constraint imposed by osmoregulation in 
A. brasiliensis, and its consequences in body size and 
maturity of females. 

Environmental sex determination (ESD) was described 
in several animal groups (West et al. 2002; Kato et al. 
2011; Gamble and Zarkower 2012). In peracarids such as 
talitrid amphipods, the effects of environmental factors on 
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Figure 5. Binomial GAMM (logit scale) of sex ratio (females/males), 
showing the partial effect of the retained environmental predictor. 
The standardized mean of the response variable (mean = 0) is 
equivalent to a prevalence of 0.57 in the original scale of 
observations. The lowest recorded value was 0.54 (in favour of 
females) at salinity = 34. The number in the y-axis label indicates the 
estimated degrees of freedom of the smooth curve. Marks on the 
x-axis indicate the measured values of each variable. 



sex ratio have been documented (Johnson et al. 2001). We 
found a prevalence of females (Fig. 5), which could be 
attributed to the local photoperiod (10-14 h of daylight). 
Studies in G. duebeni showed a female prevalence at pho- 
toperiods shorter than 14 hours (Bulnheim 1991). The 
prevalence of females decreased at higher salinities (Fig. 5). 
Three possible explanations can account for this pattern: 
(1) the infection of feminizing microsporidean protozoan 
parasites on embryos of G. duebeni populations caused a 
bias in favour of females, but its incidence is suppressed by 
salinity levels up to 25-30 (Bulnheim 1991); (2) Epigenetic 
factors like salinity could stimulate gonad maturity acting 
directly on the crustaceans' endocrine system, potentially 
determining the sexual identity of the organisms (Nagaraju 
2011); and (3) Yolk allocation could determine the sex of 
the offspring as in other zoological groups (Radder et al. 
2009). In this sense, our results found coincident trends 
for sex ratio and embryo size in relation to salinity. 
Although little is known about the mechanisms related to 
ESD (Kato et al. 2011), the models obtained for A. brasili- 
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Figure 6. Conceptual diagram of the 
connection among life-history traits in 
Atlantorchestoidea brasiliensis (see grey 
shadow). Two environmental scenarios (1 - 
high salinities, 2 - low salinities) and two 
reproductive strategies (A - low number of 
large embryos, B - high number of small 
embryos) are shown. Doted lines represent a 
shift in life-history traits for the salinity 
scenarios based on the additive nature of the 
models (see "+"). Boxes integrate the 
information between trait models for both 
salinity scenarios. The doted box indicates one 
of the possible explanations for sex 
determination (see details in text). 
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ensis suggest that larger embryos (found at low salinities) 
could produce females, while smaller embryos (found at 
high salinities) produce males (Fig. 6). 

Conclusions 

Salinity was identified as the only variable that affects all 
life-history traits of Atlantorchestoidea brasiliensis. Low- 
salinity values demand large amounts of energy for 
osmoregulation that cannot be assimilated as reserves 
(body weight) for reproductive costs. This explains why 
at low-salinity scenarios the individual weight of brood- 
ing females and the proportion of mature females were 
lower for a given size. In addition, salinity affected 
embryo size and fecundity, the latter through body 
weight. Low salinities represent a risk to offspring, being 
necessary to provide energy reserves (yolk) to each 
embryo to increase their probabilities to survive. Under 
this adverse scenario, the increase in embryo size and the 
decrease in female fecundity (due to the limited space in 
its pouch) represent a trade-off. Furthermore, yolk alloca- 
tion could explain why large embryos produce females, 
determining the high prevalence of this sex at low salini- 
ties, which should compensate low individual fecundities. 
Our results highlighted that variables acting on physiol- 
ogy influence the energy balance and determine intercon- 
nected responses of life-history traits mainly associated 
with reproduction, having consequences at the popula- 
tion level. In addition, other variables related to environ- 
mental unpredictability (in our case beach slope) can 
have significant effects on fitness through organisms' sur- 
vival. In summary, this study revealed that field data of 
multiple traits and modelling approaches that incorporate 
random factors are useful to unveil environmental drivers 
that define concerted responses of life-history traits, pro- 
viding insights into understand the potential underlying 
mechanisms. 

Acknowledgments 

We express our gratitude to the "Benthic Ecology Group" 
of UNDECIMAR for field and laboratory assistance. We 
thank Florencia Cabrera and Matthias Seaman for helpful 
grammar revisions and comments. Financial support 
provided by Fondo Clemente Estable (4034), Global Envi- 
ronmental Facility (GCP URU 030 GFF), PEDECIBA 
PhD programme and The Pew Charitable Trusts is 
acknowledged. 

Conflict of Interest 

None declared. 



References 

Barboza, F. R., J. Gomez, D. Lercari, and O. Defeo. 2012. 
Disentangling diversity patterns in sandy beaches along 
environmental gradients. PLoS ONE 7:e40468. 

Brazeiro, A., and O. Defeo. 1996. Macroinfauna zonation in 
microtidal sandy beaches: is it possible to identify patterns 
in such variable environments? Estuar. Coast. Shelf Sci. 
42:523-536. 

Bulnheim, H. 1991. Zur okophysiologie, sexualitat und 
populationsgenetik litoraler Gammaridea — ein iiberblick. 
Helgol. Mar. Res. 45:381-401. 

Defeo, O., and J. Gomez. 2005. Morphodynamics and habitat 
safety in sandy beaches: life-history adaptations in a 
supralittoral amphipod. Mar. Ecol. Prog. Ser. 293:143-153. 

Defeo, O., and A. McLachlan. 2005. Patterns, processes and 
regulatory mechanisms in sandy beach macrofauna: a 
multi-scale analysis. Mar. Ecol. Prog. Ser. 295:1-20. 

Emery, K. O. 1961. A simple method of measuring beach 
profiles. Limnol. Oceanogr. 6:90-93. 

Gamble, T., and D. Zarkower. 2012. Sex determination. Curr. 
Biol. 22:257-262. 

Gomez, J., and O. Defeo. 1999. Life history of the sandhopper 
Pseudorchestoidea brasiliensis (Amphipoda) in sandy beaches 
with contrasting morphodynamics. Mar. Ecol. Prog. Ser. 
182:209-220. 

Gomez, J., and O. Defeo. 2012. Predictive distribution 

modeling of the sandy-beach supralittoral amphipod 

Atiantorchestoidea brasiliensis along a macroscale estuarine 

gradient. Estuar. Coast. Shelf Sci. 98:84-93. 
Hastie, T., and R. Tibshirani. 1990. Generalized additive 

models. Chapman and Hall, London, U.K. 
Herrick, J. E., and T. L. fones. 2002. A dynamic cone 

penetrometer for measuring soil penetration resistance. Soil 

Sci. Soc. Am. J. 66:1320-1324. 
Hurley, D. E. 1968. Transition from water to land in 

amphipod crustaceans. Am. Zool. 8:327—353. 
Isaksson, C, B. C. Sheldon, and T. Uller. 2011. The challenges 

of integrating oxidative stress into life history biology. 

Bioscience 60:194-202. 
Johnson, W. S., M. Stevens, and L. Watling. 2001. 

Reproduction and development of marine peracaridans. 

Adv. Marine Biol. 39:105-260. 
Joshi, J., B. Schmid, M. C. Caldeira, P. G. Dimitrakopoulos, 

J. Good, R. Harris, et al. 2001. Local adaptation enhances 

performance of common plant species. Ecol. Lett. 

4:536-544. 

Kato, Y., K. Kobayashi, H. Watanabe, and T. Iguchi. 2011. 
Environmental sex determination in the branchiopod 
crustacean Daphnia magna: deep conservation of a doublesex 
gene in the sex-determining pathway. PLoS Genet. 7: 
el001345. 

Kirkwood, T. B. L., and S. N. Austad. 2000. Why do we age? 
Nature 408:233-238. 



2013 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



3923 



Linkages Among Life History Traits 



J. Gomez ef al. 



Lercari, D., and O. Defeo. 2006. Large-scale diversity and 
abundance trends in sandy beach macrofauna along full 
gradients of salinity and morphodynamics. Estuar. Coast. 
Shelf Sci. 68:27-35. 

Levitan, D. R. 2000. Optimal egg size in marine invertebrates: 
theory and phylogenetic analysis of the critical relationship 
between egg size and development time in echinoids. Am. 
Nat. 156:175-192. 

Marques, J. C, S. C. Goncalves, M. A. Pardal, L. Chelazzi, 
I. Colombini, M. Fallaci, et al. 2003. Comparison of Talitrus 
saltator (Amphipoda, Talitridae) biology, dynamics, and 
secondary production in Atlantic (Portugal) and 
Mediterranean (Italy and Tunisia) populations. Estuar. 
Coast. Shelf Sci. 58:127-148. 

McLachlan, A., and A. C. Brown. 2006. The ecology of sandy 
shores, 2nd ed. Academic Press, Burlington, MA. 

Moran, A. L., and J. S. McAlister. 2009. Egg size as a life 

history character of marine invertebrates: is it all it's cracked 
up to be? Biol. Bull. 216:226-242. 

Morritt, D., and J. I. Spicer. 1998. The physiological ecology of 
talitrid amphipods: an update. Can. J. Zool. 76:1965-1982. 

Nagaraju, G. P. 2011. Reproductive regulators in decapod 
crustaceans: an overview. J. Exp. Biol. 214:3-16. 

Normant, M., and I. Lamprecht. 2006. Does scope for growth 
change as a result of salinity stress in the amphipod 
Gammarus oceanicus 7 . J. Exp. Mar. Biol. Ecol. 334:158-163. 

Pigliucci, M. 2001. Phenotypic plasticity: beyond nature and 
nurture. Johns Hopkins Univ. Press, Baltimore, MD. 

Radder, R. S., D. A. Pike, A. E. Quinn, and R. Shine. 2009. 
Offspring sex in a lizard depends on egg size. Curr. Biol. 
19:1102-1105. 

Ramus, A. P., and R. B. Jr. Forward. 2012. The physiological 
ecology of the supratidal amphipod Talorchestia longicornis. 



Comp. Biochem. Physiol. A: Mol. Integr. Physiol. 
161:159-165. 

Ricklefs, R. E., and M. Wikelski. 2002. The physiology/ 

life-history nexus. Trends Ecol. Evol. 17:462-468. 
Schmidt-Nielsen, K. 1997. Animal physiology: adaptation 

and environment. Cambridge Univ. Press, 

Cambridge, U.K. 
Sheader, M. 1983. The reproductive biology and ecology of 

Gammarus duebeni (Crustacea: Amphipoda) in Southern 

England. J. Mar. Biol. Assoc. U.K. 63:517-540. 
Sheader, M. 1996. Factors influencing egg size in the 

gammarid amphipod Gammarus insensibilis. Mar. Biol. 

124:519-526. 

Short, A. D. 1999. Handbook of beach and shoreface 

morphodynamics. John Wiley, London, U.K. 
Spicer, J. I., P. G. Moore, and A. C. Taylor. 1987. The 

physiological ecology of land invasion by the talitridae 

(Crustacea: Amphipoda). Proceedings of the Royal Society 

of London. B. Biol. Sci. 232:95-124. 
Steele, D. H, and V. J. Steele. 1991. Effects of salinity on the 

survival, growth rate, and reproductive output of Gammarus 

lawrencianus (Crustacea, Amphipoda). Mar. Ecol. Prog. Ser. 

78:49-56. 

Truchot, J. P. 1990. Respiratory and ionic regulation in 

invertebrates exposed to both water and air. Annu. Rev. 

Physiol. 52:61-76. 
West, S. A., S. E. Reece, and B. C. Sheldon. 2002. Sex ratios. 

Heredity 88:117-124. 
Wood, S. 2011. gamm4: Generalized additive mixed models 

using mgcv and lme4. R package version 0.1-5. Available at 

http://CRAN.R-project.org/package=gamm4 



3924 



© 201 3 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



